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ABSTRACT 

^through  an  analysis  of  the  absorption  spectrum  of  Nd^  ;  Nap"  Alumina 
and  derivation  of  Judd/Ofelt  parameters,  a  model  for  the  nearest  neighbor 
crystal  field  of  the  ion  has  been  developed.  The  proposed  model  indi¬ 
cates  that  residual  sodium  in  partially  exchanged  samples  is  responsible 
for  distorting  the  symmetric  orientation  of  conduction  plane  oxygen.  This 
in  turn  induces  a  shift  in  the  Nd  ion  within  a  sixfold  symmetric  arrange¬ 
ment  of  oxygen  ions  thereby  inducing  a  linear  field  term  in  the  crystal 
field  parameters.  This  is  responsible  for  the  anomalous  absorption  strength 


INTRODUUTION 

Recent  work  has  described  how  optically  active  trivalent  transition  metal 
and  rare  earth  ions  can  be  substituted  for  monovalent  Na^  in  Sodium  Beta'' 
Alumina  (SB"A)  by  using  low  temperature  exchange  processes  <1,2>.  Jansen 
et.al.  <3>  used  this  approach  to  prepare  small  platelets  of  various  exchange 
percentages  and  observed  laser  action  at  1.06  microns  in  some  of  them.  In 
addition,  they  reported  the  existence  of  an  anomalously  strong  absorption 
line  at  5800  angstroms.  In  the  present  paper,  we  provide  a  more  detailed 
examination  of  the  absorption  spectra  by  performing  a  so-called  Judd/Ofelt 
analysis  <4,5 >  .  Judd/Ofelt  parameters  are  obtained  by  using  the  preliminary 
x-ray  diffraction  results  for  the  structures  of  sodium  and  rare  earth  beta" 
aluminas.  Based  on  this  analysis,  we  propose  a  model  for  the  nearest  neighbor 
environment  of  the  Nd^^  ion. 

ABSORPTION  SPECTRA  AND  JUDD/OFELT  ANALYSIS 
3+  4 

Nd  ions  may  be  exchanged  for  Na  in  a  full  range  of  percentages, 

typically  between  20%  of  all  Na^  exchanged  to  over  95%  Na%xchanged. 

3+  20  -3  21 

Corresponding  Nd  ion  concentrations  range  from  3.6x10  cm  to  1.8x10 

-3  20  -3 

cm  .  A  typical  sample  of  5x10  cm  concentration  yields  the  absorption 

spectrum  shown  in  Figure  1,  taken  on  a  Cary  IR-14  spectrophotometer.  There 

are  two  major  differences  between  the  spectrum  of  Nd^^  in  YAG  and  that  for 

Nd^^  in  SB' 'A.  The  first  is  that  the  individual  absorption  bands  in  Nd^^  : 

3+ 

SB"A  are  heavily  broadened  in  comparison  to  Nd  :YAG.  This  may  be  attri¬ 
buted  to  variations  in  neighbor  distances  due  to  doping-induced  strain, 
multiple  siting  or  variations  in  residual  sodium  population  due  to  defects. 
Secondly,  one  finds  that  the  absorption  band  at  the  nominal  wavelength  of 
5800  %  is  anomalously  large.  As  has  been  pointed  out  by  Judd  <4 >  ,  this 


band  in  Nd  is  strongly  influenced  by  a  linear  field  term  in  the  nearest 
neighbor  crystal  field.  In  the  current  work,  we  focus  on  the  origin  of  this 
anomalous  absorption  as  it  provides  both  the  key  to  the  understanding  of  the 
nearest  neighbor  environment  and  to  efficient  optical  pumping  at  5800 

The  Judd/Ofelt  analysis  begins  with  the  calculation  of  the  oscillator 
strength  P  for  the  transitions  between  the  initial  ^19/2  gi^ound  state 
manifold  to  the  final  states  labelled  by  RS  coupled  states  that  form  the 
major  portion  of  these  final  states.  The  oscillator  strength  for  a  group 
of  closely  spaced  transitions  is  given  by 


'’"""  +  1)1  |<i|rlf>l^ 


P  =  X  [87i‘mv/3h(2J  +  l)] 


(1) 


i.f 


where 


X  =  (n^+  2)^/9n 


is  the  correction  for  the  bulk  polarizability  for  a  medium  of  index  of  re¬ 
fraction  n,  (2J+1)  accounts  for  the  multiplicity  of  the  ground  state,  and 
<i|rlf>  are  the  matrix  elements  for  the  electric  dipole  operator  between 
initial  and  final  states.  Experimentally,  the  oscillator  strengths  are  derived 
from  the  absorption  spectra  by  the  expression  given  by  Hoogschagen  <6>  as 

p  ,  JOOOjlc! - 1  j  (2) 

where  k(X)  is  the  absorption  strength  measured  at  a  wavelength  X,  C  is  the 
molar  concentration  and  is  Avogadro's  number.  Since  the  spectrophotometer 
yields  absorption  plots  of  logd/I^)  versus  wavelength  through  a  sample  of 
thickness  I,  the  areas  under  these  curves  can  be  used  to  find  the  oscillator 
strengths  for  each  group  of  transitions  by  using  the  expression 


P  =  lP-00  !"C  X  /  log(I/I-)dX 

ire^N  \^Cl  J  ° 

0  0 

In  Table  1  we  tabulate  the  values  for  oscillator  strength  for  the  various 
transitions  for  Nd^''^  in  chlorides  <4;»YA6  <7  >  and  SB' 'A  hosts.  The  transi 
tion  at  5800  8  is  seen  to  be  an  order  of  magnitude  stronger  in  SB"A  than 
in  either  the  chloride  or  YA6  environment. 

The  group  theoretical  approach  of  Judd  and  Ofelt,  based  on  the  work  of 
Racah  <8>,  has  become  the  popular  technique  for  deriving  eigenvectors  and 
transition  strengths  for  f*^  electron  configurations  of  the  rare  earths, 
particularly  for  Nd^^  in  a  variety  of  hosts.  In  this  formulation,  the 
oscillator  strength  is  given  by 

P  =  •2rvT^^^l<i|lU^^^  11  f  >  1^ 

A 

where  <illu^^^llf>  are  the  reduced  matrix  elements  for  irreducible  tensor 
operators  between  the  initial  and  final  states  and  the  parameters  T^^^  are 
given  by 

T^^^  =  X[8/m/3h(2J+l)](2  +l)2(2t+l)'^  h  ^(X.t)  ^  |  A^pl  ^ 

The  dependence  on  the  crystal  field  environment  enters  via  the  crystal 
field  parameters  A^p  given  by 

A^p  =  (-l)^'^’i^eq^r^  Y^p[47T/2t+l ]** 

where  the  r^  are  the  individual  nearest  neighbor  positions,  the  Y^p  are 
the  sperical  harmonics  and  the  summation  is  over  the  indivdual  nearest 
neighbor  ions  with  the  Nd  ion  taken  at  the  origin.  The  variable  =  (X,t) 
contains  the  free  ion  radial  integrals  and  the  various  3-J  and  6-J  symbols. 


By  obtaining  values  for  the  oscillator  strengths  from  the  absorption  spec¬ 
trum,  a  least  squares  fit  enables  determination  of  the  T  parameters.  Again, 

3+ 

in  Table  2,  we  compare  the  values  obtained  for  Nd  in  the  three  previously 
mentioned  crystal  hosts.  Additionally,  we  note  that  Krupke<9>  has  found  it 
convenient  to  use  a  parameter  to  characterize  many  of  the  glass  hosts 
used  at  Lawrence  Livermore  Labs.  The  conversion  between  T^^^  and  is 

given  by 

T^^^  =  X[8TrV3h(2J+l)]«(j^j  (7) 

ORIGIN  OF  THE  T  PARAMETERS 

Conceivably,  given  accurate  data  on  the  positions  of  the  nearest  neighbor 

3-f 

ions  and  accurate  calculations  of  the  radial  integrals  for  the  Nd  ion,  the 
T  parameters  and  therefore,  the  oscillator  strengths  should  be  calculable 
from  first  principles.  Both  Judd  <4 >  and  Krupke  <9>  have  performed  such  a 
calculation  on  Nd^"*^  in  a  variety  of  environments  and  have  had  encouraging  but 
not  complete  success.  The  inability  to  attain  better  than  an  approximate 
match  between  experimental  and  theoretical  calculations  is  usually  blamed  on 
an  inaccurate  knowledge  of  the  radial  integrals,  usually  known  only  for  the 
free  ion.  In  the  current  work,  we  take  a  slightly  different  approach  to  the 
problem.  Preliminary  x-ray  data  <10>  on  totally  exchanged  samples  indicates 
that  rare  earth  ions  may  substitute  into  the  SB"A  system  in  sites  of  inversion 
symmetry.  This  implies  that  electric  dipole  transitions  are  forbidden  from 
parity  arguments  alone.  However,  all  of  the  samples  studied  in  the  current 
work  have  some  percentage  of  sodium  remaining  in  the  conduction  plane.  Further 
this  x-ray  work  has  indicated  that  near  neighbor  ions  to  the  dopant  can  be 
distorted  from  their  rest  position  under  the  influence  of  the  trivalent  dopant. 


Based  on  the  fact  that  such  near  neighbor  distortions  are  expected  in  the  con¬ 
duction  plane,  we  have  developed  a  model  that  treats  conduction  plane  ions  as 
parameters  that  can  be  adjusted  to  yield  a  best  fit  to  the  experimentally  ob¬ 
served  T  parameters.  Indeed,  because  of  the  matching  difficulty  found  in 
previous  work,  it  would  probably  be  best  to  match  ratios  of  the  T  parameters 
rather  than  their  absolute  values. 

One  should  always  be  skeptical  of  such  an  approach  since,  with  enough 
parameters,  a  good  theoretical/experimental  fit  can  be  obtained  with  a  possibly 
nonsensical  group  of  parameter  values.  However,  in  the  current  work,  we  are 
aided  by  the  abnormally  high  value  of  T2  and  further,  by  the  fact  that  the 
linear  A^p  crystal  field  terms  appear  only  in  T2.  It  will  be  shown  that,  to 
achieve  a  close  fit  to  the  experimentally  determined  T  parameters,  it  will 
only  be  necessary  to  assume  that  the  neodymium  ion  and  nearest  neighbor  oxygen 
ion  are  slightly  displaced  from  their  x-ray  determined  positions  in  completely 
exchanged  and  unexchanged  samples  respectively.  We  propose  that  these  dis¬ 
placements  are  simply  the  response  to  the  mutual  electrostatic  interaction 
felt  between  remaining  sodium  and  near  neighbor  oxygen. 

STRUCTURE  REVIEW  AND  CALCULATION  OF  T  PARAMETERS 

We  first  review  the  structure  of  the  SB"A  system.  As  shown  in  Figure  2, 
the  structure  consists  of  spinel  blocks  of  closely  packed  A1  and  0  separated  by 
loosely  packed  slabs  containing  Na  and  0  .  The  open  space  in  these  slabs 

permits  rapid  Na*  migration  in  two  dimensions.  Figure  3  shows  the  results  of 
projecting  the  conduction  slab  ions  onto  the  mid-plane,  or  so-called  conduction 
plane.  In  a  typical  composition  of  Na^  Mg  g^  AI^q  33  0^^,  where  Mg  sub¬ 
stitutes  for  A1  in  the  spinel  block  to  stabilize  the  beta''  phase,  sodium  ions 


occupy  5/6  of  two  crystallogrphically  equivalent  sites  (the  Beevers-Ross  and 
anti-Beevers-Ross  sites),  these  positions  being  alternately  shifted  slightly 
up  and  down  from  the  mid-plane.  A  third  crystallographic  position  (mid¬ 
oxygen  site)  is  unoccupied.  Trivalent  rare  earths  may  be  substituted  for 
the  conduction  plane  sodium  by  immersion  of  the  crystal  in  rare  earth  molten 
salts  <2>  ,  leaving  the  spinel  blocks  unchanged.  Recent  work  <10>  has  shown 
that,  when  gadolinium  is  totally  exchanged  for  sodium, gadolinium  preferentially 
occupies  mid-oxygen  sites  as  shown  in  Figure  4.  A  small  percentage  of  cells 
are  found  to  contain  gadolinium  at  the  Beevers-Ross  sites.  Unpublished  pre¬ 
liminary  work  <11  >  indicates  that  neodymium  also  occupies  the  smae  sites 
as  gadolinium  with  a  possibly  slightly  higher  percentage  occupation  of  the 
Beevers-Ross  sites. 

As  will  be  shown,  variations  in  the  positions  of  conduction  plane  0 
3+ 

and  Nd  are  crucial  to  the  Judd/Ofelt  analysis.  Fourier  synthesis  maps 

indicate  that  not  only  are  the  mid-plane  0  distorted  towards  the  mid-oxygen 

3+ 

position  by  0.5  A  but  the  Gd  position  is  found  to  vary  in  an  elliptical 
fashion  over  several  tenths  of  angstroms  <10>  . 

3+ 

Consider  first  the  occupation  of  Nd  at  the  mid-oxygen  si'*’®,  along 
with  the  nearest  neighbor  oxygen  ions  in  the  spinel  blocks  above  and  below 
the  conduction  plane  (see  Figure  5).  It  can  be  seen  that  the  Nd  occupies 
a  site  of  inversion  symmetry,  which  must  be  discounted  as  spectroscopically 
inconsistent.  However,  the  situation  changes  slightly  in  the  presence  of 
sodium.  As  the  exchange  process  is  begun,  replacing  three  sodium  ions  with 
each  neodymium  ion,  sodium  ions  are  not  expected  to  remain  as  nearest  neigh¬ 
bors  to  neodymium.  However,  as  the  exchange  process  continues,  a  situation 


; 


Is  quickly  reached  In  which  oxygen  In  the  conduction  plane  sees  an  unbalanced 
distribution  of  sodium  or  other  neodymium.  As  a  result,  conduction  plane 
oxygen  is  expected  to  shift  to  an  asymmetrical  arrangement  about  the  neodymium 
ion.  This  asymmetry  leads  to  a  high  value  for  the  linear  field  term  at  the 
neodymium  ion  and  a  correspondingly  large  value  for  T^.  However,  this  effect 
alone  cannot  account  for  the  values  of  the  observed  T  parameters;  a  shift  of 
the  Nd^^  ion  must  also  be  included.  As  the  exchange  process  is  continued  and 
sodium  is  totally  exchanged,  distortions  in  the  conduction  plane  oxygen  sites 
are  now  caused  by  adjacent  neodymium,  but  will  not  be  as  severe  as  under  the 
presence  of  sodium.  Therefore,  one  expects  the  oscillator  strength  (average 
absorption  cross  section  per  Nd  ion)  to  start  out  at  low  values  for  weak  ex¬ 
change,  to  rise  to  a  peak  at  some  optimum  exchange  percentage  and  to  decrease 
again  to  low  values  as  total  exchange  is  approached.  As  shown  in  Figure  6, 
this  is  exactly  the  effect  observed.  Indeed,  preliminary  work  indicates  that 
mobil  sodium  is  responsible  for  a  reduction  of  the  oscillator  strength  at 
increased  temperature.  Samples  that  are  nearly  totally  exchanged  show  little 
temperature  dependence  of  oscillator  strength.  However,  the  residual  oscillator 
strength  at  virtually  100%  sodium  exchange  indicates  that  sodium  alone  cannot 
be  responsible  for  the  loss  of  inversion  symmetry,  and  may  possibly  be  caused  by 
adjacent  neodymium  distorting  oxygen  positions  in  the  conduction  plane.  In  our 
parameter  fitting  procedure,  we  will  not  attempt  to  model  the  site  positions 
at  high  exchange  percentage,  but  only  the  oxygen  and  neodymium  positions 
required  to  yield  a  good  fit  with  the  observed  T  parameters  at  moderate  sodium 
exchange  percentages. 

With  this  introductory  framework  in  place,  we  now  present  the  theoretically 
calculated  T  parameters  for  the  proposed  model.  We  use  the  technique  described 


by  Judd  <4 >  and  Krupke  <9>  and  include  the  contributions  due  to  5g  configur¬ 
ations.  As  noted  above,  the  undistorted  nearest  neighbor  environment  shows 
inversion  symmetry,  therefore  T  parameters  will  be  calculated  for  various 
plausible  distortions  from  inversion  symmetry.  In  Figure  7,  we  plot  theore¬ 
tically  calculated  T  parameters  for  the  introduction  of  monovalent  sodium  at 

3+  2- 

a  given  distance  from  the  Nd  ion  in  an  otherwise  symmetric  arrangement  of  0  . 

Figure  0  displays  calculated  T  parameters  for  an  arrangement  in  which  one  of 

2 

the  conduction  plane  0  ”  ions  is  moved  a  given  distance  from  a  position  of 

3+ 

inversion  symmetry  and  finally.  Figure  9  displays  T  parameters  for  the  Nd 
ion  initially  in  a  position  of  inversion  symmetry,  after  which  it  is  displaced 
in  the  conduction  pla'i»'  along  a  line  connecting  the  conduction  plane  oxygen. 

From  an  examination  of  Figure  7,  it  is  seen  that  a  single  sodium  ion  cannot  be 
solely  responsible  for  the  experimentally  observed  T  parameters.  Table  3 
shows  the  experimentally  determined  T  parameters  for  the  same  sample  shown  in 
Figure  1.  The  second  set  of  T  parameters  in  Table  3  is  for  a  single  sodium 
ion  at  a  distance  of  4.29  angstroms  from  the  Nd  ion,  the  nearest  site  expected 
to  be  occupied  by  sodium  after  exchanging  a  single  Nd  into  the  conduction  plane. 
As  seen  from  the  table,  both  the  absolute  values  of  the  T  parameters  and  their 
ratios  are  not  at  all  a  good  match  to  the  experimental  values.  By  allowing  both 
an  asymmetrical  shift  in  conduction  plane  oxygen  (Figure  8)  and  a  shift  in  the 
neodymium  position  (Figure  9),  a  good  match  can  be  achieved  as  shown  in  the 
last  column  of  Table  3  and  in  Figure  10.  We  stress  that  the  distortions  des¬ 
cribed  in  Figure  10  are  certainly  not  uniquely  determined  by  the  T  parameters, 
rather  we  consider  the  orientation  of  ions  in  Figure  10  to  be  an  "average" 
of  many  possible  positions.  However,  we  do  suggest  that  the  asymmetric  arrange¬ 
ment  of  conduction  plane  oxygen  and  shift  of  the  neodymium  ion  are  the  type 
of  distortion  that  must  exist  to  account  for  the  observed  T  parameters. 


We  finally  consider  the  possible  contribution  due  to  Nd  at  Beevers-Ross 
sites,  shown  by  x-ray  data  to  be  occupied  a  small  percentage  of  the  time.  Table 
4  shows  the  T  parameters  calculated  for  this  situation  with  no  oxygen  or  neo¬ 
dymium  distortion.  When  the  T  parameters  are  calculated  for  this  site,  they 
must  be  multiplied  by  the  fractional  occupation  of  this  site  to  properly  account 
for  their  effect  on  the  absorption  spectrum.  When  this  is  done,  the  additional 
contriubtion  to  the  T  parameters  for  occupation  of  the  Beevers-Ross  site  can 
probably  be  ignored. 

It  must  be  emphasi2ed  that  the  success  of  the  fitting  procedure  depends 
on  the  ability  to  "tweak"  the  Na^  ion  positions  and  that  the  values  of  the  T 
parameters  are  rather  sensitive  to  these  positions.  Such  an  artificial  luxury 
could  not  be  enjoyed  by  Judd  or  Krupke  since  their  work  attempted  to  calculated 
T  parameters  for  presumably  well-known  ion  positions.  However,  we  contend  that 
the  close  fit  between  experimental  and  theoretical  T  parameter  ratios  is  strong 
support  for  the  technique  employed  in  this  work. 

CONCLUSION 

3+ 

By  analyzing  the  absorption  spectrum  of  Nd  :SB"A,  we  have  derived  the 

Judd/Ofelt  parameters  that  characterize  the  local  crystal  field  environment 

for  the  Nd^^  ion.  We  have  proposed  a  model  for  the  local  crystal  field  that 

includes  the  effect  of  residual  sodium  in  the  conduction  plane.  We  have  shown 

that  the  experimentally  determined  T  parameters  can  be  accounted  for  by  a 

3+ 

slight  displacement  of  the  Nd  ion  away  from  its  usual  mid-oxygen  site  and  by 
a  slight  displacement  of  the  conduction  plane  oxygen  away  from  a  symmetric 
arrangement  around  the  neodymium  ion. 
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TABLE  2:  T  PARAMETERS  -  From  Least  Squares  Fit 

-20 

(units  of  10“  sec.) 


NdCK:6HoO 


.125 

.87 

8.1 

1.69 

1.73 

2.2 

3.13 

3.53 

1.4 

TABLE  3.  T  PARAMETERS  -  Experimental  ve.  Theoretical 
(units  of  10"^®  eec.) 

Nd:SB'  'A  |  Single  Na'*’  0  4.3  A  I  0^”  and  Hd^  Distortion 
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Figure  1.  Typical  absorption  spectrun  for  Nd^^ :  SB* 'A,  5x10*®  /cm' 
concentration,  .26  mm  thickness. 

Figure  2.  Unit  cell  for  unexchanged  SB* 'A. 

Figure  3.  Conduction  plane  for  unexchanged  SB* 'A,  indicating  one 
of  six  possible  sodium  ion  vacancies. 

Figure  4.  Conduction  plane  for  Gd**:SB*'A,  indicating  occupation 
of  mid-oxygen  site  and  shifted  oxygen  ion  position.  Site  labels 
are  from  Reference  10. 

Figure  5.  Nearest  neighbor  environment  for  Gd^^  :SB''A  following 
the  x-ray  analysis  presented  in  Reference  10. 

Figure  6.  Experimental  variation  of  T  parameters  versus  Nd  ion 
concentration . 

Figure  7.  Variation  of  T  parameters  for  Nd*^  versus  monovalent 
ion  position  for  an  otherwise  symmetric  crystal  field. 

Figure  8.  Variation  of  T  parameters  versus  0*'  position  with 
opposite  0  held  at  2.55  A  from  mid-oxygen  site. 

Figure  9.  Variation  of  T  parameters  versus  Nd^^  shift  from  a 
position  of  inversion  symmetry. 

Figure  10.  Possible  site  distortion  for  Nd*  :SB''A.  The 
distortion  is  a  composite  of  the  shifts  described  in  Figures  8 


CONDUCTION  PLANE 
SODIUM  ALUMINA  -  UNEXCH 


Na 


unit  cell 


NTRATION  (cm'3) 


DL/413/83/01 

GEN/413-2 


TECHNICAL  REPORT  DISTRIBUTION  LIST.  GEN 


No. 

Copies 

Office  of  Naval  Research  2 

Attn:  Code  413 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 

Dr.  Bernard  Douda  1 

Naval  Weapons  Support  Center 
Code  5042 

Crane,  Indiana  47522 

Commander,  Naval  Air  Systems  1 

Command 

Attn:  Code  310C  (H.  Rosenwasser) 
Washington,  O.C.  20360 

Naval  Civil  Engineering  Laboratory  1 

Attn:  Or.  R.  W.  Drisko 

Port  Hueneme,  California  93401 


Defense  Technical  Information  Center  12 
Building  5,  Caneron  Station 
Alexandria,  Virginia  22314 


DTNSRDC  1 

Attn:  Dr.  G.  Bosmajian 
Applied  Chemistry  Division 
Annapolis,  Maryland  21401 

Dr.  William  Tolies 

Superintendent  1 

Chemistry  Division,  Code  6100 
Naval  Research  Laboratory 
Washington,  D.C.  20375 


No. 

Copies 

Dr.  David  Young  1 

Code  334 

NORDA 

NSTL,  Mississippi  39529 

Naval  Weapons  Cefiter  1 

Attn:  Dr.  Ron  Atkins 
Chemistry  Division 
China  Lake,  California  93555 

Scientific  Advisor  1 

Commandant  of  the  Marine  Corps 
Code  RD-1 

Washington,  D.C.  20380 

U.S.  Army  Research  Office  1 

Attn:  CRO-AA-IP 
P.O.  Box  12211 

Research  Tn angle  Park,  NC  27709 

Mr.  John  Boyle  1 

Materials  Branch 
Naval  Ship  Engineering  Center 
Philadelphia,  Pennsylvania  19112 

Naval  Ocean  Systems  Center  1 

Attn:  Dr.  S.  Yamamoto 
Marine  Sciences  Division 
San  Diego,  California  91232 


^7 


